Optimal Control of Brownian Inventory Models 
with Convex Inventory Cost: 
Discounted Cost Case* 



J. G. Dai^ and Dacheng Yao^ 
October 29, 2011 

Abstract 

We consider an inventory system in which inventory level fluctuates as a Brownian 
motion in the absence of control. The inventory continuously accumulates cost at a 
rate that is a general convex function of the inventory level, which can be negative 
when there is a backlog. At any time, the inventory level can be adjusted by a positive 
or negative amount, which incurs a fixed positive cost and a proportional cost. The 
challenge is to find an adjustment policy that balances the inventory cost and adjust- 
ment cost to minimize the expected total discounted cost. We provide a tutorial on 
using a three-step lower-bound approach to solving the optimal control problem \mdcr 
a discounted cost criterion. In addition, wo prove that a four-parameter control band 
policy is optimal among all feasible policies. A key step is the constructive proof of the 
existence of a unique solution to the free boundary problem. The proof leads naturally 
to an algorithm to compute the ioiiv parameters of the optimal control band policy. 
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1 Introduction 



Dai and Yao [6] studied the optimal control of Brownian inventory models under the long- 
run average cost criterion. This paper is a companion of [6]. It studies the same Brownian 
inventory models, but under the discounted cost criterion. Its main purpose is to provide 
a tutorial on the powerful, lower-bound approach to proving the optimality of a control 
band policy among all feasible policies. The tutorial is rigorous and, except the standard 
Ito formula, self contained. In addition, this paper contributes to the literature by proving 
the existence of a "smooth" solution to the free boundary problem with a general convex 
holding cost function. As a consequence, a four-parameter optimal control band policy is 
shown to be optimal. Our existence proof also leads naturally to an algorithm to compute 
the optimal control band parameters. 

The introduction in [6] gives detailed descriptions of the Brownian inventory models, 
control band policies, and the lower-bound approach. It also gives an extensive literature 
review. Most of the development there including the motivation to study non-linear holding 
cost function and the literature review applies to this paper as well and it will not be 
repeated here. In the rest of this introduction, we highlight the development that is specific 
to the discounted cost case. 

As in [6], inventory position is assumed to be adjustable, either upward or downward. 
All adjustments are realized immediately without any leadtime delay. Each upward adjust- 
ment with amount > incurs a cost K + k^, where K >Q and A; > are the fixed cost 
and the variable cost, respectively, for each upward adjustment. Similarly, each downward 
adjustment with amount ^ incurs a cost of L + £^ with fixed cost L > and variable cost 
£ > 0. In addition, we assume that the holding cost function /i : M — t- M+ is a general 
convex function that satisfies some minimal assumptions in Assumption 1. The objective 
is to find some control policy that balances the inventory cost and the adjustment cost 
so that, starting from any initial inventory level x, the (infinite-horizon) expected total 
discounted cost is minimized. When both upward and downward fixed costs are positive, 
the model is an impulse control problem. When both fixed cost are zero, the corresponding 
Brownian control problem is a singular control or instantaneous control problem. It was 
demonstrated in Section 6 of [6] that a singular control problem is much easier to solve 
than an impulse control problem and a two-parameter control band policy is optimal. This 
control band policy can be considered as the limit of a sequence of four-parameter control 
band policies, each of which is optimal for an impulse control problem. Therefore, in this 
paper, we restrict ourselves to impulse control problems; namely, we assume that K > {) 
and L > 0. Although in this paper we do not consider the singular control problem or 
the Brownian control problem when the inventory backlog is not allowed, our proof for 
the existence of an optimal control band policy for the impulse control problem can be 
extended to cover these two cases. These extensions were carried out in Sections 6 and 7 
of [6] in the average cost setting. 



1 



When the inventory holding cost function is hnear, namely, 



h{x) 



hix if X > 0, 
—pix if X < 



(1.1) 



for some constants pi > and hi > 0, Constantinides and Richard [5] proved that a 
four-parameter control band policy is optimal under the condition that 



As explained in [5] hi/f3 is the present value of the holding cost of keeping one unit of 
inventory now to infinity. If /ii//3 < i, it will never be optimal to reduce the inventory level 
as long as L > 0. Similarly, if < k, it will never be optimal to increase the inventory 
as long as K > 0. Thus, condition (1.2) is also necessary for a four-parameter control band 
policy to be optimal. Baccarin [1] sketched a proof that a four-parameter control band 
policy is also optimal when the holding cost is quadratic given by 



where hi > 0, pi > 0, /12 > and p2 > 0. In his proof, condition (1.2) is not needed 
any more as long as /12 > and p2 > 0. Baccarin [1] deferred the detailed proof for the 
existence of a solution to the four-parameter free boundary problem to an online supple- 
ment. Unfortunately, this document can no longer be located over the Internet. Assuming 
K = L > and k = i = 0, Plehn-Dujowich [10] proved that a three-parameter control 
band policy is optimal when the holding cost function h satisfies 



Both the linear cost in (1.1) and the quadratic cost in (1.3) do not satisfy the smoothness 
condition in (1.4). 

In this paper, when the holding cost function is assumed to be general, satisfying 
Assumption 1 in Section 2, we prove that a four-parameter control band policy is optimal. 
Assumption 1 on the convex holding cost function h is considerable weaker than those 
in literature. The cost functions in [1, 5, 10] all satisfy Assumption 1. Condition (2.6) 
in Assumption 1 is analogs to (1.2) and is automatically satisfied for h in (1.3). Similar 
to the companion paper [6], we adopt the three-step lower-bound approach in our proof. 
In the first step, we prove that if there exists a "smooth" test function / that satisfies a 
set of differential inequalities, the function / dominates the value function at every initial 
inventory level x. In the second step, given a control band policy, it is shown that the 



hi- I3k>0 and pi - I3£ > 0. 



(1.2) 




(1.3) 



h and h' are continuous; 

h is strictly concave and single-peaked; 

\h\, \h'\ and \h"\ are bounded by a polynomial. 



(1.4) 
(1.5) 
(1.6) 



2 



value function within the band can be obtained as the unique solution to a second order 
differential equation. In the third step, a solution to a free boundary problem is shown to 
exist and satisfy the conditions for / in the first step. 

The result in step 1 is known as the "verification theorem" in literature. All three prior 
papers [1, 5, 10] invoked the verification theorem in Richard [12], which in turn generalized 
the pioneering work of Bensoussan and Lions [2, 3]. This tutorial advocates the lower 
bound approach that was also adopted by Harrison et. al [8] and Harrison and Taksar [9]. 
The advantage of this approach is that, except for applying the standard Ito formula, it 
is self-contained, and therefore this approach can readily be rigorously adopted in other 
related settings. 

The free boundary problem is specified using the well known "smooth-pasting" method 
(see, e.g., [4]). Solving the free boundary problem in Step 3 is the most difficult task. 
We prove the existence of a solution to the free boundary problem that has four free 
parameters. Though our proof is similar to the one in [5], where a linear holding cost 
function is used, our proof is considerably more difficult. Unlike the proof in [5], our 
proof is also constructive so that it leads naturally an algorithm to compute the four 
parameters of the optimal control band. Recently, Feng and Muthuraman [7] developed an 
algorithm to compute the parameters of an optimal control band policy for the discounted 
Brownian control problem. They illustrated the convergence of their algorithm through 
some numerical examples. However, the convergence of their algorithm was not established. 

The rest of this paper is organized as follows. In Section 2, we define our Brownian 
control problem. In Section 3 we present a version of Ito formula that does not require 
the test function / be function. A lower bound for all feasible policies is established in 
Section 4. Section 5 shows that under a control band policy, the value function within the 
band can be obtained as a solution to a second order ordinary differential equation (ODE). 
Under the assumption that a free-boundary problem has a unique solution that has desired 
regularity properties. Section 6 proves that there is a control band policy whose discounted 
cost achieves the lower bound. Thus, the control band policy is optimal among all feasible 
policies. Section 7 is a lengthy one that devotes to the construction of the solution to the 
free-boundary problem. In the section, the parameters for the optimal control band policy 
are characterized. Section 7 constitutes the main technical contribution of this paper. 

2 Impulse Brownian Control Models 

Let X = {X{t),t > 0} be a Brownian motion with drift fi and variance cr^, starting from 
X. Then, X has the following representation 

X{t) = x + iJ,t + aW{t), t > 0, 

where W = {W{t),t > 0} is a standard Brownian motion that has drift 0, variance 1, 
starting from 0. We assume W is defined on some filtered probability space {Q, {Tt},T, P) 



3 



and W is an {J^t}-martingale. Thus, W is also known as an {J'^tj-standard Brownian 
motion. We use X to model the netput process of the firm. For each t >0, X(t) represents 
the inventory level at time t if no control has been exercised by time t. The netput process 
will be controlled and the actual inventory level at time t, after controls has been exercised, 
is denoted by Z{t). The controlled process is denoted by Z = {Z{t),t > 0}. With a slight 
abuse of terminology, we call Z{t) the inventory level at time t, although when Z{t) < 0, 
is the backorder level at time t. 
Controls are dictated by a policy. A policy is a pair of stochastic processes {Yi,Y2) 
that satisfies the following three properties: (a) for each sample path a; G J7, Yi{uj, •) G D, 
where B is the set of functions on M_|_ = [0, oo) that are right continuous on [0, oo) and have 
left limits in (0, oo), (b) for each u, Yi{uj, •) is a nondecreasing function, (c) Yi is adapted 
to the filtration {J-t}, namely, Yi(t) is J^j-measurable for each t > 0. We call Yi{t) and 
Y2{t) the cumulative upward and downward adjustment, respectively, of the inventory in 
[0, t]. Under a given policy {Yi,Y2), the inventory level at time t is given by 

Z{t) = X{t) + Yi{t)-Y2{t) = x + aW{t)+fit + Yi{t)-Y2it), t>0. (2.1) 

Therefore, Z is a semimartingale, namely, a martingale aW plus a process that is of 
bounded variation. 

Because K is assumed to be positive, we restrict upward controls that have a finitely 
many upward adjustment in a finite interval. This is equivalent to requiring Yi to be 
piecewise constant function on each sample path. Under such an upward control, the 
upward adjustment times can be listed as a discrete sequence {Ti{n) : n > 0}, where the 
nth upward adjustment time can be defined recursively via 

Ti{n) = mf{t > Ti(n - 1) : AYi{t) > 0}, 

where, by convention, Ti(0) = and AYi(t) = Yi{t) — Yi{t—). The amount of the nth 
upward adjustment is denoted by 

ei(n) = yi(ri(n)) - Y,{Ti{n)-) n = 0, 1, . . . . 

It is clear that specifying such a upward adjustment policy Yi = {Yi(t),t > 0} is equivalent 
to specifying a sequence of {(ri(n), ^i(n)) : n > 0}. In particular, given the sequence, one 
has 

iVi(t) 

Yiit) = YI (2-2) 

1=0 

and Ni{t) = max{n > : Ti(n) < t} is the number of upward controls [0,t]. Thus, it is 
sufficient to specify the sequence {(Ti(n), ^i(n)) : n > 0} to describe an upward adjustment 
policy. Similarly, since L > 0, it is sufficient to specify the sequence {{T2{n) , ^2{'n)) : n > 0} 
to describe a downward adjustment policy and 

N2{t) 

Y2{t) = J2 e2»- (2.3) 

1=0 
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Merging these two sequences, we have the sequence {(r„,^„),n > 0}, where T„ is the nth 
adjustment time of the inventory and is the amount of adjustment at time T„. When 
> 0, the nth adjustment is an upward adjustment and when ^„ < 0, the nth adjustment 
is a downward adjustment. The pohcy {Yi,Y2) is adapted if r„ is an {J^t}-stopping time 
and each adjustment ^„ must be ^Tn- measurable, In general, we allow an upward or 
downward adjustment at time t = 0. By convention, we set Z(0— ) = x and call Z(0— ) the 
initial inventory level. By (2.1), 

z(o) = x + yi(o)-y2(o), 

which can be different from the initial inventory level Z(0— ). 

Under a feasible policy (p = {(yi(i), l2(i)} with initial inventory level Z(0— ) = x and 
a discount rate /3 > 0, the expected total discounted cost DC(x, if) is defined to be 

DC(x,(^)=E^. / e-^^h{Z{t))dt (2.4) 

/•oo 

+ / e-^\KdNi{t)+LdN2{t) + kdYi{t) + ldY2{t)) 
Jo 

where E^, is the expectation operator conditioning on the initial inventory level being 
Z(0— ) = X. Because of (2.2) and (2.3), this Brownian inventory control model is called 
the impulse Brownian control model. Clearly, we need to restrict our feasible policies to 
satisfy 

oo 

Ex.[5]e-^^"(l + |en|)] <oo. (2.5) 

n=0 

Otherwise, T)C{x,ip) = oo. We assume the inventory cost function /i : M — t- M+ satisfies 
the following assumption. 

Assumption 1. Assume that the cost function h satisfies the following conditions: (a) 
it is continuous and convex; (b) there exists an a such that h S C^(M) except at a and 
h{a) = 0; (c) h'{x) < for x < a and h'{x) > for x > o; (d) 

lim h'{x) > ep and lim h'{x) < -k^; (2.6) 

xfoo x^—ca 

(e) h"{x) has smaller order than c^^^ as x f oo, that is 

, h"(x) , , 

lim = 0, (2.7) 

xtoo e^i^ ^ ' 



hoo 

e'^^yK\y)dy < oo, (2.5 
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where Ai = (^^ _^ 2/3cr2)i/2 _ ^ /cr2 > q. 

(f) h"{x) has smaller order than e~^^ I — oo, that is 

h"{x) 



lim 



0, 



e^''yh"{y)dy < oo. 



(2.9) 
(2.10) 



where Ag 



(^2 + 2/3^2)1/2 + ^ /^2 > 



Remark, (a) If /i is given by (1.1), (2.6) becomes (1.2), which is consistent with (13) in [5]. 
(b) The continuous and convex holding cost function h can be relaxed to be continuously 
differentiable once and twice at all but a finitely many points, (c) When limj;_s.oo h'{x) < i/3, 
it follows the same reasoning as in [5] that it will never be optimal to reduce the inventory 
level as long as L > 0. Similarly, when lim2^__j._oo h'(x) > kp, it will never be optimal to 
increase the inventory level as long as -ftT > 0. 

The following elementary lemma on the holding cost function is useful in later devel- 
opment. 



Lemma 2.1. (a) Under Assumption 1, 



li^ rxe-'^'~y-^)h"{y)dy _^ 



lim 



j^e^^^y-'''^h"{y)dy 

p{\\+\2){x-a) 



0. 



(h) Under Assumption 1, 



lim 



X2jy^^(y--^h'iy)dy 



lim 

x4-— oo 



gA2(x-a) 

Xi!^e''^^(y'^^h'{y)dy 



lim h'{x), 

x^oo 

= lim h'{x) 



g— Ai(a;— a) 

Proof, (a) We prove (2.11). The proof of (2.12) is similar and is omitted. If 



lim 

x^—oo 



(2.11) clearly holds. Now assume that 



lim 

a;4,— oo 



~>^^iy-^)h"{y)dy < oo. 



-Uy-'^)h"{y)dy = oo. 



(2.11) 
(2.12) 



(2.13) 
(2.14) 
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By using the L ' Hdpital rule, one has 



hm — 7T — —^7 ^ = nm 



xl-oo e-('^i+^2){a;-a) a;.4._oo -(Ai + A2)e-(-^i+^2)(3;-a) 

= lin, - , 

xi-oo (Al + A2)e-^2(x-a) 

= 0, 

where the last equahty is due to (2.9). 

(b) We prove (2.13). The proof of (2.13) is similar and is omitted. 

The first part of (2.6) implies that there exist a constant ci > and x" G (a, oo) such 
that for any x > x" , 

h'{x) > ci, 

which yields that 

lim / e^^^y-''^h'{y)dy > lim / e^^^^-"^^ h' {y)dy 

> ci • lim / e^2(j/-a)^y 
= oo, (2.15) 

where the first inequality is due to the assumption h'{x) > for x > a. By using the L' 
Hdpital rule, one has 

X^r e^^^y-'')h'{y)dy ^. A2e^2(^-a)^/(^) 
lim = lim 



lim h'{x). 
xtoo 



□ 



3 It 6 Formula 

In this section, we present the Ito formula, tailored to the discounted setting. 

Lemma 3.1. Assume that f E C"'^(M) and f is absolutely continuous such that f'{b) — 
f'{a) = f"{u)du for any a < b with f" locally in L^(]R). Then 

e-P'f{Z{t)) = f{Zm+ f e-P^{Tf{Z{s))-Pf{Z{s)))ds (3.1) 

rt 







+ a I e-^'nZ{s))dW{s)+ e-^^Af{Z{s)), 

"^0 0<s<t 
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where ^ 

r/(x) = -(T^/"(x) + fi.f'{x) for each j; G R such that f"{x) exists, (3.2) 

is the generator of the (/U, a'^)-Brownian motion X , and Jq e~^^f'{Z{s))dW{s) is interpreted 
as the ltd integral. 

Proof. Using (3.2) in [6] and the integration by parts formula for semimartingales (see, for 
example, Page 83 of [11]), we have (3.1). □ 

4 Lower Bound 

In this section, we state and prove a theorem that establishes a lower bomid for the optimal 
expected total discounted cost. This theorem is closely related to the "verification theorem" 
in literature. Its proof is self contained, using the Ito formula in Section 3. 
Define 

r K + k( if e > 0, 

m = { ife = o, (4.1) 

( L-l^ if e < 0. 

Theorem 4.1. Suppose that f E C"'^(M) and f is absolutely continuous with f" locally in 

L^(M). Suppose that there exists a constant M > such that \ f'{x)\ < M for all x £ M. 
Assume further that 

Tf{x) - I3f{x) + h{x) > for almost all x €R, (4.2) 

f{y)-f{x)<K + k{x-y)fory<x, (4.3) 

f{y) - f{x) <L + l{y-x) for x < y. (4.4) 

Then DC{x, (p) > f{x) for each feasible policy ip and each initial state Z{0—) = x G M. 
Proof. By Ito formula (3.1), 

rt 
'0 



-^'f{Z{t)) = f{Z{0-))+ [ e-^^{Tf{Z{s))-f3f{Z{s)))ds 

Jo 

+ a fe-^'f'{Z{s))dW{s)+ e-^'Af{Z{s)) 



0<s<t 

> /(Z(0-))- f e-^'h{Z{s))ds + a f e-^' f {Z{s))dW{s) 
Jo Jo 

+ e-^^Af{Z{s)), (4.5) 

0<s<t 
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where the inequaUty is due to (4.2). By (2.2) and (2.3), the control {{Yi{t),Y2{t)),t > 0} 
is equivalent to specifying a sequence {(T„,^„) : n = 0, 1, . . .}. Conditions (4.3) and (4.4) 
imply that and A/(Z(T„)) > — 0(^n), where (j) is given by (4.1). Therefore, (4.5) leads to 

e-^'f{Z{t)) (4.6) 

.t i-t m 

>/(Z(0-))-/ e-P'h{Z{s))ds + a e-^' nZ{s))dW{s) -Y^e-^^'-ct^iin) 
Jo Jo 



for each t >0. Fix an a; G M. We assume that 

Nit) 



/ e-^'hiZis))ds + V e-^^"<A(Cn) < oo 

■^0 n=0 y 



for each t > 0. Otherwise, T>C{x,ip) = oo and DC{x,ip) > f{x) is trivially satisfied. 
Because |/'(x)| < M, one has E^f^ e'^' f{Z{s))dW{s) = 0. Meanwhile 

f{Z{t)) < (/(Z(t))) + 

and E^[e-l^^{f{Z{t))y] is well defined, though it can be oo, where, for a 6 G M, 6^ = 
max(6, 0). Taking E^; on the both sides of (4.6) and noting /(Z(0— )) = /(x), we have 

/ /■* \ 

E,[e-^'{f{Z{t)))+]>f{x)-Ej / e-''^h{Z{s))ds + y^e-^^-^iCn)). 

V-^O n=0 ^ 



Taking limit as t — )• oo, one has 

iV{t) 



lim inf 

t— >oo 



eJ/ e-''^/i(Z(s))ds + Ve-/'^"(/>(^„)l +E,[e-^'{fiZit))y 
\Jo ^ 



>/(x). (4.7) 



The boundedness of /' implies that 

(/(x))+ <M(l + |x|), 

which further implies that 

N{t) 

{fiZ{t))y < M(l + \Zm < Mil + \x\ + + a\Wit)\ + ^iCnl). (4.8) 

n=0 

The following arguments follow the ones on Page 842 of [7]. Let z/(t) = X^^gl'^"!- Then 
(2.5) implies 



E, 



oo 



e 



-^'di^it) 



< oo. 
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From (7.5) of Taksar [13], we have 



- /3 



< oo. 



Applying Fubini's theorem, we have 

/•CO 

/ ^^E,, [ly (t)] dt < oo, 
Jo 

which, together with Lemma 4.1 of [7], imphes 

nminfe-^*E3,[i/(t)] = 0. 

Therefore, (4.8) imphes that 

hminfE^re-^*(/(Z(t)))+l< hminfE^ [e"^* (M(1 + Ixl + \fi\t + a\W{t)\ + u{t))]= 0. 



t—>oo 



□ 



5 Control Band Policies 

We use {d, D, U, u} to denote the control band policy associated with parameters d, D, U, 
and u with d < D < U < u. Let us fix a control band policy ip = {d, D, U, u} and an initial 
inventory level Z(0— ) = x. The adjustment amount ^„ of the control band policy is given 
by 

D — X if X < d, 
Co = < if d < X < u, 

U — X if X >u, 



and for n = 1, 2, 

^ ^\D-d if Z{Tn-) = d, 
\U — u if Z{Tn—) = u, 

where again Z{t—) denotes the left limit at time t, Tq = and 



Tn = inf{t > Tn-i : Z{t) € {d,u}} 

is the nth adjustment time. (By convention, we assume Z is right continuous having left 
limits.) Our first task is to obtain an expression for the value function V, where V{x) is 
the expected total discounted cost when the initial inventory level is x. We first present 
the following lemma. 
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Theorem 5.1. Assume that we fix a control band policy (p = {d,D, U,u}. If there exists 
a twice continuously differentiable function V : [d, n] — t- M that satisfies 



rV{x) - (3V{x) + h{x) = d<x<u, 



with boundary conditions 



V{d) - V{D) = K + k{D- d), 
V{u)-V{U) = L + l{u-U), 



(5.1) 



(5.2) 
(5.3) 



then for each starting point x € M, the expected total discounted cost DC{x, ip) is given by 



V{D) + K + k{D-x) for X e i-oo,d], 
V{x) = { V{x) for X e {d, u), 

V{U) + L-e{U -x) forxe[u,oo), 



(5.4) 



where V{x) is in (5.1). 

Remark. (5.2) and (5.3) imply that V is continuous at d and u. 

Proof. Consider the control band policy f = {d, D,U,u}. Let y be a twice continu- 
ously differentiable function on [d,u] that satisfies (5.1)-(5.3). Because d < Z{t) < u, by 
Lemma 3.1, we have 

r r* 1 r"^^*^ 

E,[e-'^V(Z(t))] =E,[y(Z(0))]+E,.[ e"^^ (ry(Z(s))-/3y(Z(s)))dsJ [ ^ e'^^-^ 



n=l 



where On = V{Z{Tn)) — F(Z(T„— )). Boundary conditions (5.2) and (5.3) imply that 
0n = V{Z{Tn)) - V{Z{Tn-)) = -^Cn) for u > I. Therefore, 



E,[e-^'V{Zm-^.[ViZm] 



E, 



-E, 



-E, 



^^{TV{Z{s))-(3V{Z{s)))ds 

N{t) 



Nit) 
+ E.[^ 



e-^^-On 



n=l 



[ e-^'h{Z{s))ds 
Jo 

t 

e-^'h{Z{s))ds 



E, 



E. 



n=l 

N{t) 



n=0 



Letting t — )• oo, we have 



DC(x,^) =E,.[y(Z(0))] +E,[0(eo)] 



(5.5) 
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because 

lim Ex.[e-^*y(Z(t))] = 0. 
If Z{0—) = X £ {d,u), we have Z{0) = Z{0—) = x and .^o = 0; then 

DC(x,(^) = V{x). 

\i X < d, under control band policy ip = {d, D,U,u}, Z immediately jumps up to D. 
Therefore, Z{0) = D and = D — x, then 

E,[y(Z(0))] = V{D), E.i^iCo)] = <t){D-x) = K + k{D - x), 

which, together with (5.5), implies that 

DC(x, if) = V{D) + K + k{D- x). 

The analysis for the case x > u is analogous and is omitted. □ 

We end this section by explicitly finding a solution V to (5.1)-(5.3). 
Proposition 1. Let (p = {d,D, U,u} be a control band policy with 

d < D < U < u. 

Define 

V{x) = Aie^i^ + Bie-^^'^ + Vo{x), 

where 

2 1 



\l + A2 



0-2 Al + A2 



(5.6) 



. _b2{Vo{D)-VQ{d)+K + k{D-d))-bi{Vo{U)-Vo{u) + L + i{u-U)) 
1 — r r 1 \P-') 

ai02 — 0201 

^ a2{Vo{D) - Void) + K + k{D - d)) - ai{Vo{U) - Voiu) + L + i{u - U)) . 

Bi = . (5.8) 

0261 - ai02 

Then V is a solution to (5.1)-(5.3). In (5.7) and (5.8), we set 

ai = 6^1'=' - e^i^, 02 = e^i" - e^^^, (5.9) 
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Proof. Let 

Ai = Ufi^ + 2/3^2)1/2 _ J ^ 

A2= [(^2 + 2/3^2)1/2 + 

so that z = Xi and z = — A2 are two solutions of the quadratic equation 



1 



fj^^z^ + /xz - /3 = 0. 

The homogenous ordinary differential equation (ODE) 

rg-f]g = 

has two independent solutions gi{x) and g2ix), where 



and g2ix) = e 



-X2X 



Let 



and 



= det ( ^,^5^1 = -(A, + A2)e(^-^^)- / 



ai{x) 
a2{x) 



^1 2 
- —r^g2{y)^h{y)dy = 

a W{y) (T^ 

1 2 



1 



Ai + A2 ct2 
1 2 
Ai + A2 



h{y)dy, 

Hy)dy, 



la w{y)^ Ai 

where a is the minimum point of the convex inventory cost function h. Then the non- 
homogenous ODE (5.1) has a particular solution Vb(x) given by 

2 1 



^0(2;) = [ai{x)gi{x) + a2ix)g2ix)] 



cr^ Ai + A2 



h{y)dy 



^\x(x-y) 



h{y)dy 



A general solution V{x) to (5.1) is given by 

V{x) = ^le^i^ + + Vq{x). 

Boundary conditions (5.2) and (5.3) become 

(^le^i"' + 5ie-^2rf + _ i^Aie^'^ + ^le"^^^ + yo(^)) = K + k{D -d), (5.11) 

(^,e^i" + 5ie-^2" + Vb(n)) - {Aie^^^ + Sie'^^^ + yo(f/)) = L + £{u - U). (5.12) 

Using the coefficients defined in (5.9)-(5.10), we see the boundary conditions (5.11) and 
(5.12) become 

Aiai + Bihi + Void) - V(i{D) = K + k{D- d), 
Aia2 + Bib2 + Vo{u) - Vo{U) = L + i{u - U), 

from which we have unique solution for Ai and Bi given in (5.7) and (5.8). □ 
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6 Optimal Policy and Optimal Parameters 

Theorem 4.1 suggests the following strategy to obtain an optimal policy. We hope that 
a control band policy is optimal. Therefore, the first task is to find an optimal policy 
among all control band policies. We denote this optimal control band policy by ^p* = 
{d* , D* ,U* ,u*} with the expected total discounted cost 

( V{D*)+K + k{D* -x) for X G (-oo,d*], 
V{x) = < vix) for X G ((i*,n*), (6.1) 

[ V{U*)+L-i{U* -x) for X G [u*,oo), 

for any starting point x G M. We hope that V can be used as the function / in Theorem 
4.1. To find the corresponding / that satisfies all the conditions of Theorem 4.1, we provide 
the conditions that should be imposed on the optimal parameters 

V'{D*) = -k, V'iU*)=e, (6.2) 
V'{d*) = -k, V'{u*) = L (6.3) 

See Section 5.2 of [6] for an intuitive derivation of these conditions. Under condition (6.3), 
V IS aC^ function on M. Therefore, (6.3) is also known as the "smooth-pasting" condition. 

In this section, we will first prove in Theorem 6.1 the existence of parameters d* , D* , U* 
and u* such that the value function V, defined on corresponding the control band 

policy if* = {d*,D*,U*,u*} satisfies (5.1)-(5.3) and (6.2)-(6.3). As part of the solution, 
we are to find the boundary points d* , D* , U* and u* from equations (5.1)-(5.3) and (6.2)- 
(6.3). These equations define a free boundary problem. We then prove in Theorem 6.2 that 
the function V in (6.1) with parameters d* , D* , U* and u* satisfies all the conditions in 
Theorem 4.1; therefore, the control band policy (p* is optimal among all feasible policies. 

To facilitate the presentation of Theorem 6.1, we first find a general solution without 
worrying about boundary conditions (5.2) and (5.3). Proposition 1 shows that V is given 
in the form 

V{x) = Aie^^"" + Bie-^^"" + Vo{x) for x G M, (6.4) 

where Vb is given in (5.6). Since Ai and Bi are yet to be determined, which need d* , D* , 
U* and u* , V is also yet to be determined. Differentiating both sides of (5.1) with respect 
to X, we have that 

9{x) 



is a solution to 

Tg{x)- f3g{x) + h'{x) = for ah x G M \ {a}. (6.6) 
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V'{x) 

' ' ^A. 



(6.5) 



(t2 Ai + A2 



^^(^-y)h{y)dy + Xi / e^^^^-y^h{y)dy 



g{x) in (6.5) can be rewritten as 



V'ix) 






2 1 


A2 Te-^^Cx 

a 


Ct2 Ai + A2 




2 1 




ct2 Ai + A2 



0-2 Al + A2 



A2 Ja 



A2(x-a) 



(6.7) 



where the third equahty uses the assumption that h[a) = and in the last equahty A and 
B satisfy 



1 1 Aia. ^ . 2 11 



Cr2 Al + A2 A 



Al + A2 As 



gAaa^ = -A2S1. 



(6.^ 



The following theorem characterizes optimal parameters {d* , D* ,U* ,u*) and parameters 
A* and B* in (6.7) via solution g = gA,B- 

Theorem 6.1. Assume that the holding cost function h satisfies Assumption 1. There 
exist unique A* , B* , d* , D* , U* , u* with 



d* < XI < D* < U* < X2 < u* 
such that g{x) = gA,B{x) in (6.7) satisfies 



D 



g{x) + k 



g{x)-e 

g[d) = -k, 

g{D) = -k, 

9(u) = e, 

9{u) = t 



dx = —K, 
dx = L, 



(6.9) 

(6.10) 

(6.11) 

(6.12) 
(6.13) 
(6.14) 
(6.15) 



Furthermore, g has a local minimum at xi < a and a local maximum at X2 > a. The 
function g is strictly decreasing on (— oo,xi), strictly increasing on {xi,X2) and strictly 
decreasing again on (x2,oo). 
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■k 



I 








X 



Figure 1: There exist xi < X2 such that the function g decreases in (—00, xi), increases 
in (xi,X2), and deceases again in (x2,oo). Parameters d* , D* , U* and u* are determined 
by g{d*) = g{D*) = —k, g(U*) = g{u*) = i, the shaded area between U* and u* is L, and 
the shaded area between d* and D* is K. In the interval g is the derivative of the 

relative value function associated with the control band policy {d* ,D* , U* ,u*}. 

If g satisfies all conditions (6.6), (6.10)-(6.15) in Theorem 6.1, V{x) in (6.4) clearly 
satisfies all conditions (5.1)-(5.3) and (6.2)-(6.3). The proof of Theorem 6.1 is long, and 
we defer it to Section 7. 

Theorem 6.2. Assume that the holding cost function h satisfies Assumption 1. Let d* < 
D* < U* < u* , along with constants A* and B* , be the unique solution in Theorem 6.1. 
Then the control hand policy ip* = {d* , D* ,U* ,u*} is optimal among all non- anticipating 
policies. 



Proof. Let 




for X G (—00, d*] 
for X £ {d* ,u*), 
for X G [u* ,00), 



(6.16) 



and 




x) for X G (—00, d*) 

for X G {d*,u*), 
) for X G (n*, 00), 



(6.17) 



with 



(6.18) 
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where ^^^^e'^^M* = Ai^^, ^^^^e^^a^* ^ _x^bI and Vo{x) is given by (5.6). 
Therefore, 

V'{x) = gA'',B''{x) forxeM. (6.19) 

We now show that V satisfies ah the conditions in Theorem 4.1. Thus, Theorem 4.1 
shows that the expected total discounted cost under any feasible policy is at least V{x). 
Since V{x) is the expected total discounted cost under the control band policy if* with 
starting point x, V{x) is the optimal cost and the control band policy ip* is optimal among 
all feasible policies. 

First, V{x) is in C'^{{d* ,u*)). Condition (6.10) implies 

V{d*) - V{D*) = - gA*,B- {x)dx = K + k{D* - d*) (6.20) 
and (6.11) implies 

V{u*)-V{U*)= / gA*,B*{x)dx = L + l{u* -U*). 



(6.18) implies that V satisfies 

VV{x) - fiV{x) + h{x) = {), for x G [d* , u*] . 

By Theorem 5.1, defined in (6.17) must be the discounted cost under control band policy 
if . 

Now, we show that V{x) satisfies the rest of conditions in Theorem 4.1. Conditions 
(6.12) and (6.15) imply that truncated function V'{x) is continuous in M. Therefore, 
V E Ci(M). Clearly, V'{x) = -k for x e {-oo,d*] and V'{x) = £ for x e [u*,oo). Let 

M= sup \gA*,B*{x)\. 

a;£ [d* ,u*] 

We have \ V'{x)\ < M for all x G M. Because 

TV - l3V{x) + h{x) =rV - ^V{x) + h{x) = for xe[d*,u*]. 

In particular 

TV{d*) - pV{d*) + h{d*) = 

and 

TV{u*) - ^V{u*) + h{u*) = 0. 

It follows from part (a) and part (b) of Lemma 7.2 that d* < xi < a < X2 < u* , V"{d*) = 
V"{d*) = g'{d*) < and V"{u*) = V"{u*) = g'{u*) < (see Figure 1). Thus, we 
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have fiV'id*) - (3Vid*) + h{d*) > and fiV'{u*) - f3V{u*) + h{u*) > 0. _Now, for x < d* , 
TV{x)-pV{x)+h{x) = n{-k)-P{V{d*)-k{x-d*))+h{x) > ^V'{d*)-pV{d*)+h{d*) > 0. 
Similarly, for x > u*, rV{x) - ^V{x) + h{x) = ^(^) - P{V{u*) + l{x - u*)) + h{x) > 
^lV'{u*) - ^V{u*) + h{u*) > 0. 

Now we verify that V satisfies (4.3). Let x,?/ G M with y < x. Then, 

V{x)-V{y) + k{x-y) = f [g{z) + k]dz 

AxAD*)yd* 

> / [giz] + k]dz 

J{yVd*)AD* 

> [ [g{z) + k]dz 
Jd* 

= -K, 

where the first inequality follows from g{z) = —k for z < d* and g{z) = g{z) > —k for 
D* < z < u* and g{z) = i > —k for z > u* , and the second inequality follows from the 
fact that g{z) = g{z) < —k for z G [d*,D*]; see. Figure 1. Thus (4.3) is proved. 
It remains to verify that V satisfies (4.4). For x,y G M with y > x. 

V{y)-V{x)-£{y-x) = r [g{z) - £]dz 



X 



Ayhu*)\JU* 
< / [9{z)-£]dz 

J{xVU*)Au* 



u 



< / [g{z)-£]dz 

Ju* 
= L, 

proving (4.4). □ 

7 Optimal Control Band Parameters 

This section is devoted to the proof of Theorem 6.1. We separate the proof into a series of 
lemmas. 

Since h is convex, one has h'{x) < h'{y) whenever the derivatives at a; < y exist. It 
follows that Ymix^ah' {x) and lim^.^(i /i'(a;) exist. Define 

h'{a—) = limh'{x) and h'{a+) = lim/i'(x). 

x'la x^a 

We have h'{a—) < h'{a+). Recall the function g in (6.7). Using the integration by parts. 
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one has 



9{x) 



2 1 



0-2 A1+A2 

+ Ub - h'ia-) + £e'^(y-^h"{y)dy)e- 



X2{x-a) 



+ j^h'{x) for X < a 



2 1 

0-2 A1+A2 



M ^ - h'{a+) - £'e-^i(s^-'^)/i"(y)dyje^i(^-") 



(7.1) 



+ f (5-/l'(o+)-X,''e^2{l/-a)/^//(y)^y\g-A2(x-a) +^h'{x) for X > o 



It foUows that 



/(x) = < 



2 1 

(T^ A1+A2 



■(^S - h'{a-) + /"e-^2(y-a)/j//(y)^y^g-A2(x-a) 



2 1 

(t2 A1+A2 



S - h'{a+) - £'e^2fe-a)/j//(y)^yU-A2(x-a) 



for X < a, 



for X > a 



(7.2) 



and 



/'(x) = < 



2 1 
5^ A1+A2 



(a - /i'(a-) + />-^i(s'-")/i"(y)dy)Aie^i(^-») 



+ (b - h'ia-) + /;e^2(^-'^)/i"(y)dy A2 



-A2(x— a) 



2 1 

5^ A1+A2 



A - h'{a+) - e-^i(2'-")/i"(y)dy) Aie^i^^-'^) 



+ (^B- h'{a+) - 1^ e^2(?/-a);j//(y)^y^_)^2g-A2(x-a) 



for X < a, 



for X > a. 



(7.3) 



Define 



+00 



A = h'{a+)+ e-^''^y-^h"{y)dy, B = h'{a-)- e^-'^y-^^h" {y)dy 



We have the following lemma. 

Lemma 7.1. Assume that h satisfies Assumption 1, then 

h'{a-) < A, (7.4) 

B < h'{a+). (7.5) 

Proof. Assumption 1 (c) says that h'{x) < for x < a and h'{x) > for x > a, we have 

h'{a-) = lim/i'(x) < < lim/i'(x) = h'{a+). (7.6) 

x'l'a x'4-a 

If /^+°° e-^^^y-''^h"{y)dy > 0, then (7.4) clearly holds. Now assume that 



+00 



''■y~''^h"{y)dy = 0. 
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A = Bv i-kp) 





B 



















^^^^ ^^^^^^^^ ^ 








^ ^^^^ ^.^^^ A^ 
^^^^^^^ ^ 

<^ <^ — ^ 




A = 


h'{a-) 


i 
i 

A = h'{a+) A 


= ~A 



h'{a+) 



B = h'ia-) 

g,,,eUM' B)) 
B = B 



Figure 2: (a) The shaded region is the set of {A^B) that satisfies (7.7) and (7.8). The 
unique minimum xi = xi{A, B) G (— oo, a) is well defined for all {A, B) in this region, 
(b) For each A ^ [B_\/ (-/c/3), +oo) , there exists a unique B{A) G {B_■,A^ h'{a+)) such 
that |j^^-j(xi(A, 5(^4))) = —k. The curve B = B{A) is decreasing. 



Because h"{x) > and h"{x) is assumed to be continuous on (a, oo), we have h"{x) 
for x > a. Therefore, h must be linear in x > a. This fact and (2.6) imply that 



lim/i'(x) > 0, 



h'{a+) 

which, together with (7.6) yields (7.4). 
Similarly we can prove (7.5). 



Lemma 7.2. (a) For each B satisfying 

B<B< h'{a+), 

and each A satisfying 

B <A, 

gA,B{x) attains a unique minimum in (— oo,a) at xi = xi{A,B) G (— oo,a). 



□ 

(7.7) 
(7.8) 
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A = h'{a+) 
B -f I 




A = h\a-) S,,B^x.,{A, B)) = i A ^ A 



Figure 3: (a) The shaded region is the set of all -B) that satisfies (7.8) and (7.9). The 
unique maximum X2 = X2{A, B) S (a, oo) is well defined for all (A, B) in the region, 
(b) For each B € (— oo, A A £/3), there exists a unique A{B) G (i? V h'{a—), A) such that 
9A{B),Bi^2{A{B), B)) = £. The curve A = A{B) is decreasing. 

For each A satisfying 

h'{a-) <A<A, (7.9) 

and each B satisfying (7.8), gA,B{x) attains a unique maximum in (a, oo) at X2 = X2{A, B) G 
(a, oo). 

(b) For each fixed A and B satisfying (7.7)-(7.8), the local minimizer xi = xi{A, B) is the 
unique solution in (— oo,a) to 

(A-h'{a-)+ re-^i(^-")/i"(y)dy)e^i(^i-") (7.10) 
= (B-h'{a-) + re^2(y-a)/j//(y^^y\g-A2(xi-a)_ 
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Furthermore, g'^ ^(x) < for x G (—00, xi{A, B)), g'^ ^(x) > for x G {xi{A, B), a), and 

lim gA,B{x) = +00, (7.11) 

g%B{xi{A,B))>Q. (7.12) 

For each fixed A and B satisfying (7.8)-(7.9), the local maximizer X2 = X2{A, B) is the 
unique solution in (a, 00) to 

[a - h'{a+) - J^' e-^^(^-'^)/i"(y)dy)e^^(^2-'^) (7.13) 
= (b- h'{a+) - r e^2to-»)/j"(y)c^y^e-A2(^2-a)_ 

J a 

Furthermore, g'^^{x) > for x £ {a,X2{A, B)), g'^ ^{x) < for x G {x2{A, B),oo), and 

lim gA,Bix) = -co, (7.14) 

gXB{x2{A,B)) <0. (7.15) 

Remark, (a) The set of {A, B) that satisfies (7.7) and (7.8) is the shaded region in Figure 
2. The set of (^4, B) that satisfies (7.8) and (7.9) is the shaded region in Figure 3. 
(b) Note that 

fX2{A,B) 

{A-Xi e-^i(^-'^)/i'(y)dy)e^i("2(A,B)-a) 

J a 

(■X2{A,B) 

= {A- h'{a+) - / e-^i(^-'^)/i"(y)(iy)e^i(^2(^'-^)-'^) + h'{x2iA, B)) 

J a 
fX2(A,B) 

= [B- h'{a+) - / e^2(j/-a);^//(y)^y^g-A2(^2(A,B)-a) ^ h'^x2{A, B)) 



rX2{A,B) 

= {B + \2 e^^^y-''^h'{y)dy)e-^^^''^^^''^^-"'\ (7.16) 

J a 

where the first and third equahties fohow from integration by parts, and the second is 
due to the definition of X2{A,B) in (7.13). This provide an alternative characterization of 
X2{A,B) in (7.13). Similarly, xi{A,B) has an alternative characterization. 

Proof. We only prove the existence of xi and the properties of g{x) in x G (—00, a). The 
proof for the existence of X2 and the properties of g{x) in x G (o, 00) is similar, and it is 
omitted. 

In order to prove the existence of xi, we divide B G {B,h'{a+)) into two cases: B G 
(5,/i'(a-)] and^G (/i'(a-), /i'(a+)) . 
Case 1. S G {B,h'{a-)]. 
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Note that h is convex, we have h"{x) > for aU x G M except x = a. Therefore, 
e'*'^^y^"'^h"{y)dy > is decreasing in x € (—00, a). Then for fixed B G (i?, /i'(a— )] , 
there exists an x' with x' £ {—00, a] such that 

B = h'{a-)- r e^^^y-'''^h"{y)dy. (7.17) 



We are going to prove that g'{x) is strictly increasing in x G (— oo,x') and 

lim g'{x) = -00, (7.18) 

x\.—oo 

hm(/'(x)>0, (7.19) 

x-\x' 

g'{x)>0 forxG(x',a). (7.20) 

Since ^'(x) is continuous and strictly increasing in x G (— oo,x'), (7.18) and (7.19) imply 
that there exists a unique xi with xi G (—00, x') such that 

< 0, X < xi, 
/(x) { =0, X = xi, 

> 0, xi < X < x'. 

Combining this with (7.20), we have 

< 0, X < xi, 
g'ix){ =0, x = xi, 

> 0, xi < X < a, 

from which one proves the existence of xi and properties of ^(x) in (—00, a). 

It remains to prove that g'{x) is strictly increasing in x G (— oo,x'), and that (7.12), 
and (7.18)-(7.20) hold. We first prove that g'{x) is strictly increasing in x G (— c«,x'). For 
X G (—00, x'), 

B - h'{a-) + r e^2(y-«)/i"(y)dy > 

J X 

and 

A-h'{a-)+ r e-^'^y-''h"{y)dy > B-h'{a-)+ C e-^^^y-^^h" {y)dy 



> B-h'{a-)+ r e^^'^y-^^h"{y)dy 

J X 

> 0, (7.21) 
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where the first inequahty is due to (7.8). Using (7.3), we further have that for x G (—00, x') 

2 1 
Ai + A2 

+ [b- h'{a-) + re^2(j'-")/i"(y)dy)A2e-^2(^-")" 



g"{x) 



A-h 



'(a-) + re-^i(2^-'^)/i"(y)fiy)Aie^i(^-'^) 



> 0. 

This proves g'{x) is strictly increasing in (— oo,x') 
To see (7.18), it follows from (7.2) that 



(7.22) 



lim 



9'{x) 



lim 



2 1 



(^B-h'{a-) + j e^2(^-'^)/i"(y)dy) 

a , 



a;4,-oo g-A2(a;-a) x\-cx) Ai + A2 - 

+ (A-h'ia-)+ r e 

J X 

To evaluate this limit, we first have 

lim {A-h'{a-)+ re-^i'*'(^')(^-")/i"(y)dy)e(^i+^2)(^-") 
= lim re-^i^'(^)(j'-")/i"(y)d2/- e(^i+^2)(x-a) 

x^—oo 

= 0, 

where the last equality follows from (2.11). Next, 

lim {B-h'{a-)+ re^2/i'Wfe-a)/j//(y)^y\ 
= B-B. 

Because B - B> (7.23) and (7.24) imply that 



(7.23) 



(7.24) 



lim g{x) = —00. 

x\.—ao 



(7.25) 
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To see (7.19), it follows from (7.2) that 
2 1 



limff'(x) 



> 



Cj2 Ai + A2 L 



2 1 

0-2 Ai + A2 L 
2 1 



ct2 Ai + A2 
_ 2 1 

Cr2 Ai + A2 

> 0, 



{A-h'{a-)+ / e-^i(^-")/i"(y)(iy)e^i(^'-'^) 
-(S - h'{a-) + r e^^'^y-'''>h"{y)dy)e~^^^'^'-''^ 

Jx' 

{A-h'{a-)+ [ e-^i /i"(y)(iy)e^i(^'-'') 
(^ - h'{a-) + r e^2(^-")/i"(y)(iy)e^i(^'-") 



(^-5)e 



Ai (x'—a) 



(7.26) 



where the second and last equalities are due to (7.17) and the last inequality is due to 
(7.8). 

To see (7.20), for x G [x',a), (7.17) implies that B - h'{a-) + £ e^-^^y-^'^h" {y)dy < 0, 
which plus 

A-h'{a-)+ re-^i(^-'^)/i"(2/)dy > B-h'{a-)+ P e-^^'^y-''^h"{y)dy 

J X J X 

> B-h'{a-)+ r e^^^y-^h"{y)dy 

J X 

imply that 



2 1 

^Ai + A2 



{A-h'{a-)+ / e-^^(^-'^)/i"(y)dy)e^i(^-") 



> 



\B-h'{a-)+ / e^2(^-'^)/i"(y)dy)e~^2(^~'^) 

J X 

^ ^ ^(S-/l'(a-)+ re^2fe-a);^//(y)^y^(|gAi(a-a) _g-A2(x-a)^ 



0-2 Ai + A2 

> 0. 

Case 2. B £ {h' (a-) , h' {a+)) 

It is similar to prove (7.22), (7.25) and (7.26), we have 

g"{x) > for X G (—00, a), 
lim g'{x) = —00 

x^—oo 



(7.27) 
(7.28) 



and 



lim g'{x) 



2 1 

^Ai + A2 



-(A-B) > 0. 



25 



Therefore, there exists a unique xi such that 




Limit (7.11) can immediately be obtained by (7.25) and (7.28). Inequalities (7.22) and 
(7.27) and the definition of xi easily imply (7.12). □ 

Lemma 7.3. Suppose A and B satisfy (7.7)-(7.8), for fixed B, the local minimizer xi{A, B) 
is continuous and strictly decreasing in A; for fixed A, the local minimizer xi{A,B) is 
continuous and strictly increasing in B. Suppose A and B satisfy (7.8)-(7.9), for fixed B, 
the local maximizer X2{A,B) is continuous and strictly increasing in A; for fixed A, the 
local maximizer X2{A, B) is continuous and strictly decreasing in B. 
Furthermore, 

lim xi(yl,S) = -oo, (7.29) 

limxi(yl,S) = a for A<h'ia+), (7.30) 

B-\A 

limx2(vl,S) = a for B > h'{a-), (7.31) 

limx2{A,B) = oo. (7.32) 

Proof. The Implicit Function Theorem implies the continuity of Xi{A,B), i = 1,2. Apply- 
ing the Implicit Function Theorem to (7.10) and (7.13), we have that 

< 0, (7.33) 

> 0, (7.34) 



dxi{A,B) 


2 I ^Xi{xi{A,B)-a) 


dA 


cy^ A1 + A2 g"{xi{A,B)) 


dxi{A,B) 


2 1 ^~X2{xi(A,B)-a) 


dB 


ct2Ai + A2 g"{xi{A,B)) 


dx2{A,B) 


2 1 ^\i{x2(A,B)-a) 


dA 


A1 + A2 g"{x2{A,B)) 


dx2{A,B) 


2 I g-A2(x2(A,B)-a) 


dB 


a^\i + X2 g"{x2{A,B)) 



> 0, (7.35) 

< 0, (7.36) 

where in obtaining (7.33) and (7.34) we have used g"{xi{A,B)) > in (7.12), and in 
obtaining (7.35) and (7.36) we have used g" {x2{A, B)) < in (7.15). 

Fix A satisfying A < h'(a+), when xi t a, (7-10) gives that B '[ A. From the mono- 
tonicity between xi and B, we must have (7.30). 
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We next prove (7.29). From (7.10), we have 

(^A - h'{a-) + Te-^i /i"(y)dy)e(^i+^2)(xi-a) 



= B-h'{a-) + / e^'^y-''h"{y)dy. (7.37) 

We will show that 

Urn (a - h'{a-) + r e"^i(^~")/i"(y)dy)e(^i+^2)(xi-a) ^ q_ (7 gg) 

2:14.-00 V J^^ ) 

This, together with (7.37), imphes that 

= Hm (5-/i'(o-)+ r e^^^y-''h"{y)dy) 
= hm B-h'{a-)+ [ e^''^y-"'^h"{y)dy 

xii-txi J-00 

= hm B-B, 

xi4,— 00 

from which one has that B \, B when xi ^ —00. Using the monotonicity between xi and 
B (see (7.34)), we must have (7.29). 

It remains to prove (7.38). To see this, 

lim (A-h'{a-)+ re-^i(s'-'^)/i"(y)(iy)e(^i+^2)(xi-a) 



Xll- 

= hm / e-^i(^-")/i"(y)d?/ • e(^i+^2)(xi-a) 
xii-00 J^^ 

= 0, 

where the last equahty is due to (2.11). Therefore, we have proved (7.38). 

The proof for (7.31) and (7.32) is similar. □ 

Lemma 7.4. (a) For each 

Be(-oo,AM^), (7.39) 



there exists a unique 

A{B) e {BVh'{a-),A) 

such that 

gA(B),BMA{B),B))=L (7.40) 
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Furthermore, for B G (—00, A A iP), 

Therefore, function A = A{B) is strictly decreasing in B ^ (—00, A A iji); see Figure 3 for 
an illustration. For A € {A{B),A), 

gA,B{x2{A,B))> L (7.42) 

(h) For each 

Ae {BV {-kp),+oo), (7.43) 

there exists a unique 

B{A) G {B,AAh'{a+)) 

such that 

9AMA)iMA^BiA))) = -k. (7.44) 
Furthermore, for A £ {B\/ (—kf3), 00), 

^-^(^) = _^g{Ai+A2)(xi(AB(A))-a) ^ Q 

dA \i ■ ^ ■ ^ 

Therefore, function B = B{A) is strictly decreasing in A G {BV (— A;/?),oo); see Figure 2 
for an illustration. For B G {B_,B{A)), 

gA,B{xiiA,B)) <-k. (7.46) 

(c) The two curves {{A{B),B) : B G (-00, 0]} and {{A,B{A)) : A G [0,oo)} have a 
unique intersection point (A™*,i?™*) that satisfies 

B{A'''^) = B'""^ and A{B'''') = (7.47) 

with 

< ^(0) < A'""' < A{B) < A, (7.48) 

> B{{)) > S'"* > B(A) > B. (7.49) 

See Figure 4 for an illustration. 
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A = Bv i-kfi) B " 




A = h'{a-) g,,{x,{A,B)) = i A = A 

Figure 4: The two curves {{A{B),B) : B G [B,0]} and {{A,B{A)) : A G [0,^4]} have a 
unique intersection point that satisfies < ^(0) < < A{B) < A and > 

B{0) > > B{A) > B. For any {A, B) in the shaded region g^^^^j^^{xi{A,B{A))) < -k 
and gA(B)A''^iMB),B))>£. 



Proof, (a) First, fix a i? that satisfies (7.39). We consider the value oi gA,B{x2{A, B)) for 
Ae {BV h'{a-),A). 

oA ' oA o"^ Ai + A2 Ai 



1 1 



-e 



\i{x2{A,B)-a) 



0-2 Ai + A2 Ai ' 

> 0. (7.50) 



Next we will prove that 



\iTngA,B{x2{A,B))> i (7.51) 
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and 

lim gA,Bix2iAB)) <i, (7.52) 

Al{BVh'{a-)) 

from which one has that there exists unique A(B) G (i3 V h'{a—),A) such that 

gA{B)MMMB),B))=£ 

and for A G {A{B),A) 

gA,B{x2{A,B))> t 

The derivative (7.41) foUows from the Implicit Function Theorem, being applied to (7.40). 
First we prove (7.51). Then (7.32) implies that 

lim(S + A2 / e^2(^-")/i'(y)(iy)e~^2(^2(A,B)-a) 

AtA' J a 

= lim As / e^''^y-''^h'{y)dy ■ e-^2(x2(A,B)-a) 

A-tA Ja 

^ A2j>^^(^-")^^(j/)dj/ 

xtoo gA2(x-a) 

= lim h'{x), (7.53) 



where the last equlity follows from (2.13). Equalities (7.16) and (7.53) yield that 

nX2iA,B) 

_ 'A-\i I 

Therefore, using the expression in (6.7) for g, we have 



rx2{A,B) 

lim(A - Ai / e-^i(^-'')/i'(y)dy)e^i(^2(AB)-a) ^ ^iui h'{x) 

A^A Ja 2;too 



\i-mgA^B{x2{A,B)) 
A\A 

^ ^ [(f + f )lim/i'(x) 



fJ Ai + A2 Ai A2 a;too 
= — lim h' [x) 

where the second equality uses A1A2 = ^ and the last inequality is due to the first part of 
(2.6). 
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It remains to prove (7.52). Next we consider two cases: B G {h'{a—),A A if3) and 
B € (-oo,/i'(a-)]. If B £ {h'la-),AAei3), Huiaib X2iA, B) = a in (7.31) implies that 



hmgA.MMAB)) = ^^-^[(f + ^)(i?-/^'(a+))]+^/i'(a+) 
AiB cr^ Ai + A2 Ai A2 P 

B 
J' 

Because B < £(3, we have limA^B gA,B{x2{A, B)) < L 

On the other hand, if G (— 00, /i'(a— )], (7.16) imphes that 



(7.54) 



hm iB + \2 r ' e^2(s/-a);^/. w N -A2(x2(AS)-a) 

hm (^-Ai / e-^i(j'-'^)/i'(2/)dy)e^i(^2(A,B)-a) 



<0, 

where the inequahty is because h'{a—) < by (7.6). Using the expression in (6.7) for g, 
we have 

hm gA,Bix2iAB))<0<i. 
Aih'{a-) 

(b) For a fixed A that satisfies (7.43). We can prove similarly that for B € {B,A A 
h'{a+)) 

dgAMMA,B)) ^ 2 1^ 1 ^_,^(,^(^,^)_,) ^ ^ 
(95 0-2 Ai + A2 A2 



and 



lim gA,B{xiiA, B)) = ^ lim h'ix) < -k, (7.56) 

_B4,B p xj,— CO 

where the inequality is due to the second part of (2.6). 
If ^ G (5 V {-kl3),h'{a+)), we have 

hmgAMMAB)) = ^. (7.57) 

Because A > —k(3, we have liuiB^A gA,Bixi{A, B)) > —k. Then (7.55) and (7.56) imply 
that there exists a unique B{A) E {B_, A) such that 



9AMA)iMAB{A))) = -k 
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and for B G {B,B{A)) 

gA,B{xi{A,B)) < -k. 

If ^4 G [/i'(a+), oo), we have 

lim gA,B{xi{A,B))>0>-k. 

Then (7.55) imphes that there exists a unique B{A) G {13, h'{a+)) such that 

5A,iJ(A)(^l(^'^(^))) = 

and for B £ {B,B{A)) 

gA,B{xi{A,B)) < -k. 

Applying the Imphcit Function Theorem to (7.44), we also have (7.45). 

(c) First consider the curve {{A{B),B) : B G {—oo, A A £/3)} that is determined by 
equation gA{B),B{x2{A{B), B)) = i. Consider two points 

(4(0), 0) and {A{B),B) 

on the curve {{A{B),B) : B G {—oo,A A i(3)} (see Figure 4). By part (a) of this lemma, 
we have 

A{B) < A. (7.58) 

Next we show that 

A{0) > 0. (7.59) 

To see this, (7.54) implies that 

lim 5^,0 (a;2(^,0)) = < 

from which and (7.50), one has (7.59). 
Similarly, consider two points 

(0,5(0)) and {A,B{A)) 

on the curve determined by g^j^^y^-^{xi{A, B{A))) = —k. Similar to (7.58) and (7.59), by 
part (b) of this lemma, we have 

B < B{A) < B{0) < 0. (7.60) 
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Therefore, the point {A,B{A)) is on the right side of the curve gA{B),B{x2iA{B), B) = i 
and point (0, B{0)) is on the left side of the curve. The continuity and monotonicity of the 
two curves imply that there is a unique point 

at which the two curves intersect. See Figure 4 for an illustration. It is clear from Figure 4 
that (7.48) and (7.49) hold. 

□ 

Let 

G = {{A, B) : A{B) <A<A, B<B< B{A)} (7.61) 

be the shaded region in Figure 4. The region G has four corners. They are (^™*,i?™*), 
{A,B{A)), {A,B_) and {A{B),I3). Its boundary has four pieces: the top, the right, the 
bottom and the left. 

For {A, B) e G, we have 

qaA^M, B)) < -k, gA,Bix2{A, B)) > i. (7.62) 

It follows from part (b) of Lemma 7.2 and (7.62) that there exist unique d{A, B), D(A, B), 
U{A, B) and u{A, B) such that 

d{A,B) < xi{A,B) < D{A,B) < U{A,B) < X2iA,B) < u{A,B), 
gAMdiA B)) = gAMD{A, B)) = -k, 
gAAU{A,B)) = gA,B{u{A,B)) = I, 
g'^^s{d{A,B))<0, g\j,{D{A,B))>Q, 
gXBiUiA,B))>0, <7a,bKA5))<0. 

For each {A, B) G G, define 

j-D{A,B) MA,B) 

Ai{A,B)= [gA,B{x) + k]dx, and A2{A, B) = [gA,B{x) - e]dx. 

Jd{A,B) Ju{A,B) 

Although {A, B{A)) is not in G for A G (A™*, A), these points are on the upper boundary 
of G, and 

A2{A,B{A)) 

is also well defined for A G (A™*,^). 

Lemma 7.5. There exists a unique A = Ai £ {A™^,A) such that 

A2(Ai,B(A,)) = L, (7.63) 

and for Ae (^i,^), 

A2{A,B{A))>L. (7.64) 



33 



B 






A 

[A{0),0) 


1 

int 


4i 

^A{B) 


> 


(0, ^(0) ) 

^int 
^1 ► 

B(A) 

\ j 1 / 










A 

^AM,B{A,)) = L 

g,,{x^{A,B)) = -k 








(a{b) , b)\ 

g,_,(x,{A,B)) = I 


B = B 

A = A 



Figure 5: The point {Ai,Bi) is the unique point on the top boundary of G such that 
A2(Ai,i?i) = L. For any point {A, B) on the top boundary and to the right of {Ai,Bi), 
K2{A,B) > L. 

Proof. When A goes to {A, BiA)) goes to {A''\ B'"^*). Then the definition of (^'"*, 5^°*) 
in Lemma 7.4 implies that 

hm U{A,B{A))= hm u{A,BiA)) = X2iA''\B'''^). 

Therefore, 

hm A2{A,B{A)) = 0. (7.65) 
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Fix A e (^'°*, A). One has 
dK2{A,B{A)) 



dA 

du{A,B{A)) 



[9amaM^^B{A)))-1] 



^ ' JU(A,B{A)) 
-MA)) dg^ j^^^^ix) 



UiA,B{A)) dA 



1 r'-^»»[i^..,_, 



a Ai + A2 Ju{A,B{A)) ^1 ^2 dA 

n -, fuiA,B{A)) -I -, _ 



O" Ai + A2 Ju{A,B{A)) ^1 ^1 

> 0, (7.66) 

where the second equahty is due to g^j^^j^^{U{A, B{A))) = g^-gf^^'j{u{A, B{A))) = i, the 

forth equahty is from (7.45), and the inequahty is due to u{A, B{A)) > U{A, B{A)) > 
xi{A,B{A)). Therefore A2{A,B{A)) is increasing in A G (A'°*,]4). 
We win show next that 

limA2{A, B{A)) = 00. (7.67) 
AfA 

It foUows from (7.65), (7.67) and the monotonicity of A2{A, B{A)) that there exists unique 
Ai e (^'°*,A) such that (7.63) and (7.64) hold. 
To prove (7.67), note that (7.45) imphes that 

dg,B(Ax2iA,B{A))) 



dA 



1 ^ ^ r ^ Ai(x2(A.B(A))-a) , 1 dB{A) ,^(^^(^^B(A))-a)l 

a2Ai + A2^Ai ^ X2 dA J 

2 1 r 1 \, t^^f A ~Rf A\\-„\ 1 



gAi(x2(A,B(A))-a) _ _^g(Ai+A2)(a;i(A,B(A))-a)g-A2(x2(A,B(A))-a)l 



0-2 Ai + A2 '-Ai " Ai 
> 0, (7.68) 

where the second equality is due to g'^j^^^^{x2{A, B{A))) = 0, and the inequality is due to 

X2(^,^(^)) >Xl(A^(^)). 
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For A G {A^'^^,A), {A,B{A)) on the right side of the curve gA,B{x2{A B)) = £ and 
therefore 

gA,BiA)i^2{A,B{m > (7.69) 

Fix an A' e (A™*,]4) and let 

Ml = [gA>^B^A')i^2{A',B{A')))-e)/2. 

It follows from (7.69) that Mi > 0. Then (7.68) implies that for each A £ (A', A), 

9A,BiA)iMAB{A)))> g^,j^^^,^{x2{A',B{A')))=e + 2Mi>i + Mi. 
Therefore, for each A £ (A' , A), there exist unique Ui{A, B{A)) and ui{A, B{A)) such that 
Ui{A,B{A)) < X2{A,B{A)) < Ui{A,B{A)) 

9AMA)(Ui{A,B{A)))=g^^^^^^{ui{A,B{A)))=£ + Mu (7.70) 
5!4,B(A)(f^i(^'^(^))) > 0' 9'A,B^A)(''^iA,B{A))) < 0. 

The properties of gA,B in Lemma 7.2 imply that for A G (A' , A), 

U{A,B{A)) < Ui{A,B{A)) < X2{A,B{A)) < ui{A,B{A)) < u{A,B{A)). 

This implies that 

\imui{A,B{A)) > \imx2{A,B{A)) 

AtA AtA 

> limx2{A, B'""^) 
A'tA 

= oo, (7.71) 

where the second inequality holds because (7.36) and the equality is due to (7.32). There- 
fore, for Ae {A',A), 

_ fuiA,B{A)) 

A2iA,BiA)) = / _ [5^5(^)(x)-£]dx 

JU{A,B{A)) 
.ui{A,B{A)) 



> 



Ui{A,BiA)) 

> Mi{ui{A,B{A))-Ui{A,B{A))). 
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Applying the Implicit Function Theorem to g^^^^-^{Ui{A, B{A))) = £ + Mi, we have that 
dUi{A,B{A)) 



dA 

2 1 
^ A1+A2 



2 1 

0-2 A1+A2 



J_pXi{Ui(A,B{A))-a) : ±dB(A) X2{Ui{A,B{A))-a) 
Ai ^ X2 dA ^ 

J_gAi(l/i(AB(A))-a) _ J_g(Ai+A2)(xi(AB(A))-a)g-A2(C/i(AB(^))-") 
Ai Ai 



a'AMA)^U,iA,BiA))) 



<o, 



where the second equality is due to (7.45), and the inequality is due to Ui{A, B{A)) > 
xi{A,B{A)) and g'^.j^^^^{Ui{A,B{A))) > 0. Thus, for any A E (A', A), 

Ui{A,B{A)) < Ui{A',B{A')). 
Therefore, for any A G {A', A), 

A2{A,B{A)) > Mi{ui{A,B{A)) - Ui{A' ,B{A'))), 
which, together with (7.71), implies (7.67). □ 
Define 

Bi = B(Ai). (7.72) 

It follows from (7.45) that 

B < B(A) <Bi< B""* < 0. 
where the last inequality is due to (7.49). See Figure 5 for the point {Ai,Bi). 
Lemma 7.6. (a) For B G (^, ^i], there exists unique A*{B) G [j4i,j4) such that 

K2{A*{B),B) = L. (7.73) 

(h) For Be {B,Bi], 

dA*{B) Af(e-^2(«(A*(B),S)-a) _ g-A2(t/(A*(B),B)-a)) 



dB ;^2(gAi(M{A*{i?),B)-a) _ gAi(C/(A*(B),B)-a)) 



< 0. (7.74) 
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Proof, (a) For B = Bi, Lemma 7.5 has showed that 

A*(Bi) = Ai. (7.75) 
For B G iB,Bi) and {A, B) € G, we first have 



dA Ju{A,B) dA 



U{A,B) r, 1 



'U{A,B) ^ Ai + A2 Ai 

> 0. (7.76) 

From the definition of A{B) in (7.40), we have 

hm UiA,B)= lim u{A,B)= lim X2iA,B) = X2{A(B),B). 
AiA{B) AiA{B) AiA(B) v_v / / 

Therefore, for a fixed B G (i?, i?i), 

lim A2(^,S) = 0<L. (7.77) 

Next for B G {B,Bi), we consider two cases depending on whether B G (S, i?(A)] or 
B £ {B{A), Bi). See Figure 6 for an illustration. 

We first assume that B G {B{A),Bi). For a fixed B G {B{A),Bi), by the monotonicity 
of B{-) in (7.45), there exists an A{B) G {Ai,A) such that (^(i?),i?) is on the upper 
boundary of G. It follows from (7.66) and the definition of Ai in Lemma 7.5 that 

A2iAiB),B) = A2{A{B),B{A{B))) > A2(Ai,B(Ai)) = L, 

which, together with (7.76) and (7.77), implies that there exists a unique 

A*{B)(^{A{B),A{B)) 

such that (7.73) holds. 

Now assume that B G {B, B{A)]. Following the proof for (7.67), one can prove similarly 
that 

lim A2(^, -B) = 00, 

A^A 

which, together with (7.76) and (7.77), implies that there exists a unique A*{B) G {A{B),A) 
such that (7.63) holds. By (7.41) and (7.75), we have for B G {B,Bi), 

A*{B) >Ai. 

(b) Applying the Implicit Function Theorem to A2{A*{B),B) = L, we have (7.74). 

□ 
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^^int^ ^int j A = Al A = A{B) 



A = A 



B = Bx 



B = B{A) 



B = B* 



B = B 




.g,,{x,UB)) = -k 



A,(A*(B*),B*) = -K 
AM*(B),B) = L 



g.AxM, B)) = 



Figure 6: For B G {B, Bi], there exists a unique A*{B) £ [^i, ^) such that A2iA*{B),B) = 
L. There is a unique B* G {B,Bi) that satisfies Ai{A* (B*), B*) = -K. 

For each B G {B,Bi), Lemma 7.6 shows that {A*{B),B) G G. Thus, 

gA*iB),BMA*{B),B))<-k 

and 

Ai{A*{B),B) = / + A;]dx 

Jd{A*{B),B) 

is weh defined. 

Lemma 7.7. There exists a unique B* with B* G {B,Bi) such that Ai{A* (B*), B*) = -K. 

Proof. We only need to show that Ai{A*{B), B) can take any value in (— oo,0) for B G 
{B, Bi) and is strictly increasing in B. 

It has been shown in Lemma 7.5 that that A*{Bi) = Ai and {Ai,Bi) is on the upper 
boundary of G (the blue curve in Figure 6). Therefore 



5Ai,Bi(^i(^1'^i)) = 



and 



lim gA*^B),B{xiiA*{B),B)) = g^,,j^.j^^{xiiA*{Bi),Bi))=g^^j^^ixiiAuBi)) = -k. 

B^Bi 



(7.78) 



39 



It follows that 



We now prove 



lirn Ai{A*{B),B) = 0. 



^mA,{A*{B),B) = -^. (7.79) 



First, we prove 

dgA'iB),B{M^*iB),B)) 

dB 

To see this, for B € {B, Bi), 

dgA*{B)M^M*{B),B)) 



> 0. (7.80) 



dB 

2 1 r 1 dA*{B) . 



{xiiA*{B),B)~a) , J_„-A2(a;i(A*(B),B)-a)l 



r " M{xi{A"{B),B)~a) , J_ 

^ ^ 1 r 1 Ai{e e J Ai{xi(A*(B),B)-a) 

Cj-^X^ + X^lX^ ;^2(^gAi{«{A-(B),i?)-a) _ gAi(C/(A*(B),B)-a)) 

A2 ' 

where the second equality follows from (7.74). Using the Lagrange Mean Value Theorem, 
there exist yi G {U {A* (B) , B) , u{A* (B) , B)) and 7/2 G (C/(^*(S), S), u(A*(fi), 5)) such 
that 

g-A2KA'(B),B)-a) _ ^-X2(U(A*iB),B)-a) ^ _ x^^- X2{yi~ a) ^u{A*{B),B) - U{A*{B),B)), 

(7.81) 

^\i{u{A*{B),B)-a) _ ^x^{u{A'{B),B)--a) ^ Xic^^'^y^" ''\u{A* {B) , B) - U {A* (B) , B)) . {7 .82) 
Therefore, for B G {B,Bi), 

dgA*(B)M^iiA*{B),B)) 
dB 

= ^ ^ r ^ e~^'^^'~°^ ^Aifa(A*(B).i?)-a) ^ ^g-A2(xi(A*(i?),B)-a)l 
0-2 Ai + A2 '- A2 e^i(2^2-a) A2 

2 -A2(yi-a)^_gAi(xi(A*(B),i?)-?;2) _^g-A2(2;i(A*(B),B)-yi)j 



0-2 Ai + A2 A2 

> 0, 



where the inequality holds because xi{A*{B),B) < D{A*{B), B) < U{A*{B), B) < yi and 
xi{A*{B),B) < D{A*{B),B) < U{A*{B),B) < y2- Thus, we have proved (7.80). 
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Fix an B2 G iB,Bi). Define 

M2 = ^ ' . 

2 

It follows from (7.78) and (7.80) that g^,^^,^-^^^^{xi{A* (B2),B2)) < -k and thus M2 > 0. 
Prom (7.80), one has that, for B G (^,^2), 

gA*(B),B{MA*iB),B)) < g^.^j^^^ j^^{xi{A*(B2),B2)) = -k - 2M2 < -k - M2. 

Therefore, for B e {B,B2), there exist unique di{A*{B),B) and Di{A*{B),B) such that 

di{A*{B),B) <xi{A*{B),B) < Di{A*{B),B), 
9AHB),B{dM*iB),B)) = gA.(^B),B{DiiA*{B),B)) = -k - M2, 
9'a*(B),b{'^i{A*{B),B)) < 0, > 0. 

The properties of gA*(B),B iii Lemma 7.2 imply that for B £ {B,B2), 

d{A*{B),B) < di{A*{B),B) < xi{A*{B),B) < Di{A*{B),B) < D{A*{B),B). 
Therefore, for B € iB,B2), 

^■D{A*iB),B) 



Ai{A*iB),B) = / [9A'iB)M^) + k]dx 

Jd{A*{B),B) 
r-Di{A*(B),B) 

< / [gA*{B),B{^) + ^]'^^ 

Jdi (A*(B).B) 



di(A*{B),B) 

< -M2{Di{A*{B),B)-di{A*{B),B)). 



By (7.29), (7.33) and A*{B) > ^'"S one has that 



lim xi{A*{B),B) < lim 3;i(^""*, B) = -00. (7.83) 

BX.B B\.B_ 



Because di{A*{B),B) < xi{A*{B),B), (7.83) implies that 



Now we prove 



lim di{A*{B),B) = -00. (7.84) 



hmDi{A*{B),B)> -00, (7.85) 

B^B 



which, together with (7.84), implies that 

^raA^{A*{B),B) 



< lim -M2{Di{A*{B),B)-di{A*{B),B)) 

B\.B_ 



-OO, 
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proving (7.79). 

To prove (7.85), noting the definitions of yi and y2, we have that for B G {B, Bi), 
dDi{A*{B),B) 



dB 

_ 2 1 
^ A1+A2 



1 dA*iB) x,(D,(A*(B],B)-a) , J_„-A2(Di(A*(B),B)-a) 
Xi dB ^ ^ A2^ 



9'a^^b),b(Di{A*{B),B)) 



1 



-X2{Di(A*{B),B)-a) 



A2 



^2 X, + X2g'^,^^^,,{D^{A*iB),B)) 

1 Af(e-^^(-(-^-(^)-^)-") -e-^^(^(-^*(^)-^)-")) 

Ai Xj(e^iHA*(B),B)-a) _ gAi{C/(A*(B),i?)-a)) 

1 e-^afa-'^) Ai(Di(A*(B).g)-a) , 1 ^~X2(D^(A*(B),B)-a) 
2. i. X2 e^i(!'2-'') ^ ^ A2^ 

^'Ai + A2 g'^,^^^,,{DM*{B),B)) 

2 1 1 w ,e^i(^i(^*(^)>^)-2^2) - e-^2(Di(yl*(B),B)-j/i) 
-A2(i/i-«) 



^2Ai + A2A2 <7:i.(B),B(^l(^*(^)'^)) 

<0, 

where the inequahty is due to Di{A*{B),B) < U{A*{B),B) < yi, Di{A*{B),B) < 
U{A*{B),B) < y2 and g'^,(^^^^^{Di{A*{B), B)) > 0. Therefore, we have proved (7.85). 
Finahy we show that 9M{A^^b),b) ^ ^ foUows from (7.74) that 

dAi{A*iB),B) 
dB 

Jd(A'{B),B) oB dB ' ' 

-^?™^[....),.(.(^*(i.),^)) + .] 
2 1 r'^'^''^'''\ldA*{B).,, , 1 



ct2 Ai + A2 



r ''"''"iLdAm^X.i.-a) ^ l,-A2(.-a)l^, 
Jd{A*{B),B) Ai dB A2 



2 1 i-D{A"{B),B) ^ ^2f^^-X2{u{A*{B),B)-a) _ ^-X2{U{A'-(B),B)-a)^ 

^ X1 + X2 J diAHB),B) ^AT A2(e^i(«(^*(^)'S)-«) -e^i(f^(^*{^).^)-'^)) ' ^ " 

+ _Lg-A2(x-a)l^^ 
A2 

2 11 1 



(j2 Ai + A2 A| e-^i("(^*(^)'^)-'*) - eMU{A'-{B),B)-a 
^^-X2iu(A*(B),B)~a) _ g-A2(C/(A*(B),B)-a)^^gAi(D{A-(B),i?)-a) _ ^Xr{d{A* {B),B)-a)^ 

_(^^-X2{D{A*{B),B)-a) _ ^-~X2{d{A'- {B),B)-a)^^^X^{u{A'- {B),B)-a) _ ^X^(U{A*{B),B)-a)^^ 
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If the expression inside the bracket is positive, we must have — ^ > 0. Note that 

d{A*{B),B) < D{A*{B),B) < U{A*{B),B) < u{A*{B),B). Thus, the positivity of the 
expression is equivalent to 

^-\2{u{A*(B),B)-a) _ ^-\2(U(A*{B),B)~a) ^-\2{D{A*{B),B)-a) _ ^-MidiA* (B),B)-a) 
gXi{u{A*{B),B)-a) _ g\i{U{A*{B),B)-a) ^ gAi{_D(A*(_B),_B)-a) _ gXi{d{A*{B),B)-a) ' ■^^) 

Using the Lagrange Mean Value Theorem, there exist zi G {d{A*{B),B),D{A*(B),B)) 
and Z2 G {d{A*{B), B), D{A*{B), B)) such that 

^-X2{D(A''{B),B)-a) _ ^~X2{d{A'iB),B)-a) ^ -X2e-^-^'^-''\D{A* {B) , B) - d{A* {B) , B)) , 
^\^{D{A*{B),B)-a) _ ^\^{d{A*{B),B)-a) ^ AiC^^^^^-a) ^) _ ^ _ 

Using (7.81) and (7.82), we have that inequahty (7.86) is equivalent to 

g-A2(yi-a) g-A2(zi-a) 
gAi(j/2-a) ^ gAi(^2-a) ' 

which is further equivalent to 

g-A2(s/i-2i) ^ gAi(y2-Z2)^ (7.87) 

Inequality (7.87) holds because yi > U{A*{B),B) > D{A*{B),B) > zi and y2 > U{A*{B),B) > 
D{A*{B),B) > Z2 imply that 

yi- zi> 0, y2- Z2> 0. 
Therefore, we have proved ^^^(^^^^)>^) > completing the proof of the lemma. □ 
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